Heavy-light decay topologies as a new strategy 
to discover a heavy gluon 



Cesare Bini", Roberto Contino'^, Natascia Vignaroli'^''' 



(N 



O 

CN " Dipartimento di Fisica, Sapienza Universita di Roma and IN FN, Italy 

^ c3 ^Department of Physics and Astronomy, Iowa State University, Ames, lA 50011 

0\ 
<N 

Oh 

I Abstract 

^ We study the collider phenomenology of the lightest Kaluza-Klein excitation of the 

'— ' gluon, G* , in theories with a warped extra dimension. We do so by means of a two-site 

^ effective lagrangian which includes only the lowest-lying spin-1 and spin- 1/2 resonances. 

^ We point out the importance of the decays of G* to one SM plus one heavy fermion, 

that were overlooked in the previous literature. It turns out that, when kinematically 
^ allowed, such heavy-light decays are powerful channels for discovering the G*. In 

*0 particular, we present a parton-level Montecarlo analysis of the final state Wtb that 

follows from the decay of G* to one SM top or bottom quark plus its heavy partner. 

We find that at ^/s = 7TeV and with 10fb~^ of integrated luminosity, the LHC can 
, I discover a KK gluon with mass in the range Mq* = (1.8 — 2.2) TeV if its coupling to 

a pair of light quarks is ga^qq = (0.2 — 0.5)g3. The same process is also competitive 
•fh for the discovery of the top and bottom partners as well. We find, for example, that 

r> the LHC at ^/s = 7 TeV can discover a 1 TeV KK bottom quark with an integrated 

luminosity of (5.3 — 0.61) fb~^ for gG*qq = (0.2 — 0.5)g3. 



1 Introduction 



One of the robust predictions of theories with a warped extra-dimension [T] and fields prop- 
agating in the bulk [2]-^ is the existence of heavier copies of the gluon with ~ TeV mass, its 
Kaluza-Klein (KK) excitations. The phenomenology of the lightest of the KK gluons, which 
we will denote as G*, has been extensively studied |5-10 . In the theoretically attractive 



framework where the hierarchy of masses and mixings of the Standard Model (SM) quarks 
follows from the geography of wave functions in the bulk, while the 5D Yukawa couplings are 



flavor anarchic [3||4, 11 , the couplings of the heavy gluon G* to the SM quarks are dictated 
by an inverse hierarchy, where the largest coupling is to the top quark. In this case the decay 
rate of G* to pairs of top quarks is typically large and dominates over the decay rate to two 
light quarks. 

Experiments at the Tevatron and the LHC have been testing this theoretical picture 



performing several searches for the heavy gluon G* in dijet flnal states 12,13 and in the tt 
channel 14 19]. In the latter case, some analyses do not attempt to reconstruct the top 



quarks, and are in fact searches for WbWb resonances. Others instead look for highly boosted 
top quarks and make use of dedicated top-tagging techniques to reduce the SM background, 
following the original suggestion of Refs. |5||6]. Taken all together, these experimental results 
impose important constraints on the value of the mass and couplings of the heavy gluon, 
and in turn on the parameter space of the extra-dimensional theories. 

In addition to the KK gluon, new heavy fermions are also naturally predicted in theories 
with a warped extra dimension as the Kaluza-Klein excitations of the SM quarks and leptons. 
In particular, in models where the Higgs doublet is a pseudo Nambu-Goldstone boson of a 



larger symmetry breaking 20 ,21 , the lightest KK top regulates the quadratically divergent 
contribution of the SM top quark to the Higgs mass term. By naturalness arguments, it is 
thus expected to have a mass below ~ 1 TeV. The coupling of G* to the heavy fermions is 
typically large, so that, if kinematically allowed, the heavy gluon decays to a pair of heavy 
fermions with a large rate. In this limit the total decay width of G* becomes very large, due 
both to the large value of the couplings and to the multiplicity of available decay channels, 
making its discovery quite challenging |22j. 

On the other hand, G* can also decay to one SM plus one heavy fermion. Speciflcally, 
the most important channels in scenarios with flavor anarchy are those with one SM top or 
bottom quark plus its heavy copy. If the threshold to two heavy fermions is kinematically 
closed, the branching ratio of these heavy-light channels can be large in a wide portion of 
the parameter space. At the same time the branching ratio to pairs of SM quarks is reduced, 
making the standard analyses of dijet and ti flnal states less constraining. The importance 



of the heavy-light decays at the Tevatron was already pointed out in Ref. 23 in the context 
of a different theoretical framework. No detailed study however has been done so far, to our 
knowledge, to assess the prospects of observing the heavy-light decays of G* at the LHC. 

In this paper we present a flrst analysis of this kind and outline a simple cut-based strategy 
to maximize the signal signiflcance over the SM background. Some preliminary results have 
already been presented in Ref. [24j. We flnd that, thanks to their distinct topology, decays 
to one SM plus one heavy fermion are extremely clean channels with a strong potential for 



^Here we list only those experimental searches that presently impose the strongest bounds, omitting 
previous or less constraining analyses. 
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the discovery not just of G*, but of the KK top and bottom as well. Our analysis suggests 
that the discovery reach on these latter can be larger than the one obtained by exploiting 
the more studied processes of QCD pair production 25,22,26 28 and electroweak single 



production 29 32 



Although our study is performed in the specific context of warped extra-dimensional 
theories, we believe that our results apply to a broader class of models beyond the SM where 
both a heavy spin-1 color-octect and heavy top or bottom fermions exist, like for example 



the top-quark seesaw model 33 , 34 



The paper is organized as follows. In Section |2] we define the effective two-site model that 
we adopt to study the phenomenology of the heavy gluon and heavy fermions. We derive the 
G* production cross section and the relevant decay rates, showing that the branching ratio 
of G* to one SM plus one heavy fermion is large when kinematically allowed. In Section [3] 
we study the prospects of observing the heavy-light decays of G* in Wtb events at the LHC. 
We perform a Montecarlo simulation of the signal and the main SM backgrounds and design 
a simple set of kinematic cuts that maximizes the discovery significance. In Section |4] we 
discuss the results obtained and draw our conclusions. Finally, Appendix |A] contains a few 
results used in discussing the effective model, while details on how the statistical errors are 
computed are given in Appendix |Bj 



2 An effective theory of the heavy gluon and the top 
partners 

Instead of considering the full set of particles and interactions arising in a specific model, 
we will work in the framework of an effective theory that includes only the lowest-lying 
resonances and reproduces, in the same spirit of chiral lagrangians in QCD, the low-energy 
limit of a large set of warped extra-dimensional theories with a custodial symmetry in the 
Specifically, we will construct the effective lagrangian by following the rules of 

and adopt the 'dual' language of strongly-interacting 
In this perspective, the new heavy particles arise 



bulk 35 



deconstruction defined in Ref. 36 



dynamics to describe the theory 37,38 



as composite states of a new strong sector, which are coupled to the SM elementary ones via 
linear mixing terms, leading to a scenario of partial compositeness 39 , 36 . The Higgs is also 



a bound state of the new dynamics and has direct couplings only to the composite fermions. 



2.1 The model 

We consider the case in which the composite sector has an SU{3)c x 0(4) x U{l)x global 
symmetry, with 0(4) D 5*0(4) ~ SU{2)l x SU{2)r, where 0(4) contains a discrete parity 
Plr under which SU{2)l and SU{2)r are exchanged. The two building blocks of the model 
are the elementary sector and the composite sector. The particle content of the elementary 
sector is that of the SM without the Higgs, and the SU{3)c x SU(2)l x U{1)y elementary 
fields gauge the corresponding global invariance of the strong dynamics, with Y = Tf^ + X. 
The composite sector comprises a heavy gluon, G*, which transforms as (8, l,l)o under 
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SU{3)c X SU{2)l X SU{2)r x U{1)x, the composite Higgs 



H 



:i,2,2i 



(1) 



and the following set of vector-like composite fermions: 



Q = (3,2,2) 



2/3 



T T, 



5/3 



B T, 



Q' = (3,2,2y_i 



/3 



2/3_ 



5^1/3 T'' 
-B-4/3 -B' 



f =(3,1,1)2/3 



S = (3,l,l)_i/3 



(2) 



Other composite states, for example spin-1 resonances with electroweak quantum numbers, 
are in general present in the spectrum of realistic models, but will not be considered here for 
simplicity. The quantum numbers of the composite fermions and the Higgs under SU{3)c x 
SU{2)l X SU{2)r X U{\)x are those specified in eqs.([T]),(|2]). The fermions, in particular, 
can be arranged in two fundamentals of 5*0 (5), and in fact our effective model describes the 
low-energy regime of the SO{b) / S0{4) theories introduced in Refs. 40,41 , in the limit in 



which only the leading terms in an expansion in powers of the Higgs field are retained (see 
Refs. |42||43], for two- and three-site effective theories where the full Higgs non-linearities 
are included). For this reason, we will dub our two-site model as TS5. As it will be even 



more evident shortly, it extends the two-site model of Ref. 27 to include the heavy gluon 



and the composite fermions needed to give mass to the bottom quark. 

The lagrangian that describes our effective theory is the following (we work in the elec- 
troweak gauge-less limit, and omit the terms involving the SU (2)^ x U{1)y elementary gauge 
fields, which play no role in our analysis): 



elementary 



composite 



mixing 



(3) 



elementary 



composite 



1 



G*G 



- Y,u Ti{QU} f - Y,D Tr{Q'H} B + h.c. 



(4) 



(5) 



1 M"^ 

C^xma = ^^{G,-Glf - Ail ql (T, B) - A^^ ql (T', B') 

^ 9*z (6) 

- /\tR inf - Am IrB + h.c. , 



3 



where ViTi'^'H) is the Higgs potential, and the derivative is covariant under SU{3)c trans- 
formations. The superscript i in eq.Q runs over the three SM famihes {i = 1,2,3), with 
q'l = {tL^bi), u|. = tji, d\ = bji- In addition to the Higgs mass and the self-couphngs 
contained in the Higgs potential, the lagrangian ^ has thirteen free parameters: four com- 
posite masses in the fermionic sector (mg, mg/, rhf, rh^); one for G* (Mg*)', four mass 
mixing terms (A^i, Al2, ^bR)', two composite Yukawa couplings {Y^:U, Kkd); and two 
gauge couplings (geis, 9*3)- 

By construction, the elementary fields couple to the composite ones only through the 
mass mixing lagrangian Cmixing- This implies that the SM Yukawa couplings arise only 
through the couplings of the Higgs to the composite fermions, and their mixings to the 
elementary fermions. We consider the case in which the strong sector is flavor anarchic |3| 
|4 , 11 . Under this assumption the hierarchy in the masses and mixings of the SM quarks 



follows from the hierarchy in the composite/elementary mixing parameters. Reproducing 
the full spectrum of quark masses and the Cabibbo-Kobayashi-Maskawa matrix requires 
introducing in the lagrangian three families of heavy fermions and three sets of mixing terms 
with coefficients {A^^^, A}^2) ^ur^ ^dj?}) for each SM flavor i. In this case, on the other 
hand, the mixing parameters of the light elementary quarks are small and their effect can 
be neglected in studying the G* phenomenology, so that one can focus on just the third 
generation of composite fermions. These two simplifications have in fact been applied in 
writing eqs.([5),([6]), and are justified a posteriori by our assumption of flavor anarchy in the 
strong sector, j The 'anarchic' scenario has been extensively studied in the framework of 5D 



warped models, see Refs. 44 -49 . Although flavor-changing neutral effects are parametrically 
suppressed by the degree of compositeness of the SM fermions, a mechanism which has been 
dubbed RS-GIM flavor protection |45|, it has been shown that numerically strong bounds 
on the masses of the KK gluon and of the KK fermions still generically follow from CP- 



violating observables in the kaon sector, namely ex 46 ,47 and e'/e 49 , the neutron electric 



dipole moment 45 , and b ^ S'j 48 . Several solutions have been proposed to evade these 
bounds, without renouncing to the explanation of the SM flavor hierarchy, in which the 
strong sector {i.e. the bulk of 5D models) is assumed to be invariant under additional flavor 
symmetries (50[]53] or to preserve CP (Toj. In these scenarios the KK gluon can have a mass 
as light as a few TeVs and its phenomenology is not qualitatively modified. In the following 
we will assume that some mechanism is at work to alleviate the flavor bounds and we will 
estimate the LHC reach on the heavy masses without imposing any restriction on them. 

We assume A/^2 ^ ^li, which can naturally follow, for example, from the RG flow in 
the full theory Eol. In this limit, up to small 0{Al2) effects, the Plr parity of the strong 



sector is not broken by the coupling of the elementary b^ to the latter, and the vertex Zbb 



^In fact, once produced, heavy fermions of the first two generations wiU also decay mostly to tops and 
bottoms, since flavor-changing transitions are not suppressed in the strong sector, while the couplings to the 
light SM quarks are extremely small, see the discussion in Ref. (SGj. In this sense, the addition of the first 
two generations of heavy fermions in the lagrangian ([3| would not qualitatively alter the phenomenology 
described in this paper. 

■^If the strong sector is not anarchic, the phenomenology can change significantly and additional signatures, 
like for example the unsuppressed decay of G* to light quarks, can appear, see for example Ref. |10 . Our 
analysis, on the other hand, can still be relevant although the decay modes considered here might not be 
the dominant ones. 
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is protected by large tree- level corrections |54]. |^ In fact, this protection was the main 
motivation to require the strong sector to be invariant under Plr and choose the fermionic 
representations of eq.([2]). We further require AtR ~ Af,/?, so that the small ratio mi,/mt 
follows from the hierarchy A12/ An ^ 1. Notice also that for the two-site model 



of Ref. 27 is contained as a subsector of the TS5 model. 

In order to diagonalize the two-site lagrangian ([s]) one has to make a field rotation from 
the elementary/composite basis to the mass eigenstate basis [36]. The only exception is 
given by the first two generations of elementary quarks, since they do not mix with the 
composite fermions and can thus be directly identified with the corresponding SM states. 
We performed the field rotation analytically before EWSB by treating as a small pa- 
rameter and including only the leading order terms in the perturbative expansion. The final 
lagrangian is shown in eqs. (|27|)-(|30| of Appendix |A| We use an economical notation where 
symbols denoting elementary (composite) fields before the rotation now indicate the SM 



(heavy) fields. More details will be reported in a separate publication 56 . Here we just list 
the relevant parameters that control the spectrum and the interactions of the mass eigen- 
states and discuss the numerical assumptions we made on them. The elementary/composite 
rotation in the fermionic sector can be parametrized in terms of six mixing angles, which 
can be conveniently defined as follows: 

, ^tR . Ali Al2 /„x 

tanv9tij = ^— , taxiipbR = — — , tany^L = ^ — , S2 = - — cosy^L, (") 
ruf niQ mqi 



plus the two mixing parameters defined in eq.(31) of Appendix |A[ Here sm.'^tR-, sinyjb/j, 
smifL respectively denote the degree of compositeness of t^, hR and {tL^hi)- In the gauge 
sector, the elementary and the composite couplings Qeiz, 9*3 determine the rotation angle 6*3 
and the SM gauge coupling g^, as follows: 

tan 63 = — , 5(3 = g^is cos 63 = g^^ sin 6*3 . (8) 

9*3 

The physical masses are given by 

rhf rhn ffin 
ruf = , mo = , rriT = tub = , rriT' = ms' — mni , (9) 

cos (ftR cos LpbR cos LfL 

"^Tg/a = rriTy^ = rriT cosy^L , "^s_4/3 = "^b_i/3 = "^t' , (10) 
where only the first four are independent, and 

Mg* = ^. (11) 
cos 63 

There are thus thirteen independent angles and physical masses (in addition to the Higgs 
mass and self-couplings), corresponding to the initial number of parameters of the la- 
grangian ([3]). 

^It has been also pointed out in Ref. that in the case of S'0(5)/50(4) theories the Plr parity is 
accidental at low energy if the fermions transform as fundamental representations of 5*0(5). This ensures 
automatic protection to Zbb even if the strong dynamics is not Plr invariant at the fundamental level. 
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After the EWSB, the SM top and bottom quarks acquire a mass, and the heavy masses in 
) get corrections of order {vY^/ffiY, where v = 246 GeV is the electroweak scale. 
In the following, we assume r = [vY^/rh) ^ 1 and compute all quantities at leading order 
m r. This is in fact a required condition in order to consistently neglect operators with a 
higher number of Higgs fields in the two-site effective lagrangian. We thus work under the 
hypothesis of parametrically small r, although in some of the benchmark points that we will 
consider r is numerically not very small. In those cases we expect that our phenomenological 
analysis still gives an accurate description at the qualitative level, although our formulas can 
receive quantitatively important corrections. At leading order in r, the masses of the SM 
top and bottom quarks read 

V V 

m = -^*v sin v^L sin v?ti? , = —Y^r)S2 sin (pbR . (12) 

Fixing rrit, rrih and to their experimental values gives three conditions on the set of 
parameters. As a further simplification, in what follows we will set 

Y^u = y*D = y* , sin iphR = sin , ttit' = rriT = = ruf , (13) 

so that there are five free parameters left: three mixing angles {sin (pi, sinyjj/j, tan ^3) and 
two masses [rrif, Mq*)- Alternatively, one can also trade two of the mixing angles, sin 
tan 6*3, for the two couplings F*, (7*3. 



eqs.g, (10 



2.2 Phenomenology of G" 



Starting from the rotated lagrangian (27), one can derive the phenomenology of the heavy 
gluon. The couplings of G* to the fermions are reported in eqs. (|32|)-(|35|) of Appendix |A| In 
particular, G* couples to the light quarks with strength (—(73 tan 6^3), and in general with a 
much stronger coupling to the top, bottom and heavy fermions. Gauge invariance forbids 
a coupling with two gluons and one G* at the level of renormalizable operators. The main 
rate of production at the LHC hence comes from the Drell-Yan process pp ^ qq ^ G*. 
Figure [l] shows the relative cross section for a reference value tan 9^ = 0.44 (corresponding 



to 5^*3 = 3). The analytic expression of the G* decay partial widths is given in eqs. (36) 
(39) of Appendix |a| When kinematically allowed, G* can decay to two SM fermions, ipip 



qq + tt + bb, where q = u,d, c, s; two heavy fermions, XX? where x denotes any of the heavy 
fermions; or one heavy plus one SM fermion, ipx + x4' = Tti + Bbi + TIr + BbR + c.c..[\ 
The plot on the left of Fig. |2] shows the relative branching ratios as functions of Mq* for 
reference values nif = ITeV and tan 6*3 = 0.44, simptR = 0.6, = 3 (hence simpi = 0.56). 
The corresponding total decay width of G* is shown in the plot on the right of the same 
figure, for the same reference values of parameters. 

When the decays to one heavy fermion are kinematically forbidden, Mq* ^ mf, the heavy 
gluon can decay only to pairs of SM quarks. The relative importance of the various channels 



^ More specifically, we compute the SM top and bottom masses and all decay widths of the heavy fermions 
at 0{r), and we neglect the O(r^) EWSB corrections to the masses of the heavy fermions. 

^In the most extreme case under consideration Y^, — 3, fh — 1 TeV, so that r ~ 3/4. 

^We classify all particles according to their SU{2)l x U{\)y quantum numbers, and compute all decay 
widths at 0{r). See also footnote [Sj 
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Figure 1: Cross section of the Drell-Yan production of G* at the LHC, pp ^ qq ^ G* , for 
tan^s = 0.44 ((7*3 = 3). For different values of tanO^ the cross section scales as (tan^a)^. 
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Figure 2: Left plot: branching ratios for the decay of G* to two SM fermions, ipip, two heavy 
fermions, xXj &iid one SM plus one heavy fermion, ipx + X^i & function of Mq*. Right plot: 
total decay width of G* as a function of Mq* ■ Both plots are done setting the other parameters to 
the following reference values: nif = 1 TeV and tan^s = 0.44, simptR = 0.6, = 3. 
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is controlled by tan 6*3 and the top degrees of compositeness sin (ptR, sin cpL. For small values 
of tan 6*3, which are naturally implied by the hierarchy of couplings geis ^ g*3, or large top 
degree of compositeness, the dominant channel is tt. For example, for a fully composite t^, 
simptu = 1, and tan^^3 = 0.2, Y^, = 3 (hence simpi = 0.33) one has BR{G* tt) = 0.98, 
BR{G* — i- qq) = 0.012. In this limit, ti is the best channel to discover the G*, despite the 
huge QCD background [5|[^. On the other hand, large branching ratios to pairs of light 
quarks can be obtained even for moderate top degrees of compositeness if tan 6*3 is not too 
small. This is a consequence of the cancellation that can occur between the two terms in 



eq.(33) due to the relative minus sign. For example, for sinLptR = 0.6, tan6'3 = 0.44, K,, = 3 
one has BR{G* ti) = 0.18, BR{G* — qq) = 0.69. In this case the strongest discovery 
reach (or exclusion power) comes from the dijet searches. 

For large values of Mg*, the heavy gluon mostly decays to pairs of heavy fermions, due 
to the large number of available channels and the large couplings. This is clearly illustrated 
by the plot on the left of Fig. [2J In this case a first threshold to T5/3T5/3 + T2/3T2/3 opens 
up at Mg* = 2mTg^g = 2mfCosipL — 1.66 TeV, and a second one at Mg* = 2mf = 2TeV. 
In this limit the large value of the total decay width that follows for the large multiplicity 
of decay channels (see the right plot of Fig. |2]) makes a discovery of G* very challenging. 

There is however a third scenario, which is realized for mf < Mg* < 2mT^^^. In this 
case the threshold for decaying to two heavy fermions is kinematically closed, so that the G* 
total width remains small, but decays to one SM fermion plus one heavy fermion are allowed 
and have a large branching ratio, see Fig. |2} Such heavy-light decays are possible only if the 
SM fermion and the heavy fermion have the same quantum numbers under SU {2)lxU (1) y, 
that is, the heavy fermion must be one of the partners of the top or bottom: x = T, 5, T, B. 
Due to their large Yukawa coupling, these heavy states mostly decay to one top or bottom 
plus one longitudinally-polarized vector boson or Higgs boson. The analytic expression of 



the corresponding decay widths is reported in eqs.(40)-(42) of Appendix A In the limit 



3> Tnt,mw, the branching ratios of the possible decay modes are: 
BR (T ^ ZLtR) ^BR{T ^ Mr) ~ 50% , 
BR {B WitR) ~ 100% , 

(14) 

BR{T WLhi) ^ 50% , BR{T ^ ZLti) ^ BR{T ^ Ml) ~ 25% , 

BR{B ZLbR) ^BR{B ^ Mr) ~ 50% . 

Considering the pattern of decays, we can identify five final channels: 

Bh, ft Wtb , Tt, ft Ztt, Tt, ft Mi, Bb Zbb , Bb Mb . (15) 

Their relative importance is shown in Fig. [3] as a function of Mg* , for the same reference 
values of parameters chosen for Fig. |2] One can see that each of the heavy-light topologies has 
a sizable branching ratio in a large range of G* masses above nif and below 2'mT^^^ , while the 
ti channel is suppressed compared to its value at Mg* < ruf. Notice that although BR{G* — t- 
qq) gets suppressed as well, it is still sizable for the benchmark values of parameters chosen 
in Fig. [3j While large branching fractions for the heavy-light decays are mainly implied by 
the kinematics and as such are a robust prediction, the value of BR{G* — )■ qq) is strongly 
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Figure 3: Branching ratio of the various final states that follow from the decay of G* , as functions 
of Mg* ■ The other parameters are set to the same reference values used for Fig. [2j rrij, = 1 TeV 
and tan 63 = 0.44, sin iftji = 0.6, = 3. The decay channel htt (hbh) has the same branching ratio 
of Ztt (Zbb). 




Figure 4: Branching ratios of G* to the various final states as functions of AIq*- Compared 
to Fig. |3j two different set of parameters have been chosen: simptR = 1, tan^s = 0.2, y<< = 3, 
simpL = 0.33 (left plot), and simptR = 0.8, tan ^3 = 0.2, = 3, sinLpi = 0.41 (right plot). Both 
plots are done fixing nif = 1 TeV. 
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dependent on tan 9^ and can thus be easily made small. This is clearly illustrated in Fig. |4| 
where the branching ratios to the various channels are shown for two different choices of 
parameters. The left plot refers to the same benchmark point adopted in Ref. js], sin^j^^ = 1, 
tan ^3 = 0.2, Ki, = 3 {simpL = 0.33); in this case the tt channel largely dominates over the 
others until the threshold Mq* = 2mT^^^, as a consequence of the full degree of compositeness 
of the right-handed top. The right plot refers instead to a similar set of parameters where, 
however, the degree of compositeness of is slightly smaller: simptR = 0.8, tan ^3 = 0.2, 
K = 3 {sin ip I = 0.41). It is evident how in this case the tt branching ratio is substantially 
reduced in the range ruf < Mq* < 2mT^^^, and those of the heavy-light channels, especially 
Wtb, are sizable. Notice that in this case the branching ratio of qq is extremely small, as 
due to the small value of tan6'3. 

Given their numerical importance for rrif < Mq* < 2mT^^^, the heavy-light decay chan- 
nels represent a new promising strategy for discovering the G*, in a limit in which the latter 
is still a relatively narrow resonance. The reason is mainly twofold: on the one hand, one can 
use the peculiar topology of the heavy-light decays, where the invariant mass of a subsystem 
reconstructs the physical mass of the heavy fermion, to efficiently reduce the QCD back- 
ground; on the other hand, the tt channel has a much smaller rate for rrif < Mq* < 2mT^^^, 
hence it is much less efficient as a discovery channel. Quite remarkably, as it will emerge 
from our analysis, the heavy-light decays seem extremely competitive channels for discover- 
ing the heavy fermions as well. As a final bonus, they also give the opportunity to measure 
some important features of the fermionic sector, like for example the couplings of the heavy 
fermions to the SM vector bosons, and thus determine its origin. 

Ideally, one could look for the presence of (fermionic) resonances in all the channels of 
eq.(15). Simple considerations on the size of the SM backgrounds suggest that the most 
promising final states could be the following (we assume that the Higgs boson is light and 
mostly decays to bb): Wtb; Zti,hti — )■ bbti; Zbb,hbb — )■ 46. A heavy fermion with charge 
-1-2/3 will thus appear as a (bbt) resonance in bbtt events, and possibly as a (Wb) resonance 
in Wtb events if it is a singlet of SU{2)l (T). The non-observation of this latter signal would 
in turn point in favor of a classification of the heavy top as part of a doublet of SU{2)l (T). 
A heavy fermion with charge —1/3, on the other hand, will appear as a (36) resonance in 
46 events only if it is an electroweak singlet {B) and the SM b^ is not too elementary (this 
was our assumption when we derived the estimates in eq.(14)). The (46) channel can thus 
give a way to test the degree of compositeness of 6/j. If 6r is very much elementary, on the 
other hand, even a singlet B will decay mostly to Wtb events, similarly to a doublet B. In 
this case one will observe a {Wt) resonance in Wtb events. 

The above discussion shows how a detailed picture on the quantum numbers of the 
fermionic resonances can be obtained by exploiting all the decay final states. In what follows 
we will consider the Wtb channel, which seems to be the most promising one, and we will 
study it in detail. 



3 Analysis of the Wth channel at the LHC 

In this section we discuss the prospects of observing the signal in the Wtb channel at the LHC. 
We will present a simple parton-level analysis aimed at assessing the LHC discovery reach. 



10 




and focus on final states with one lepton. Our event selection will be cut-based and driven 
by simplicity as much as possible. We consider two center-of-mass energies: y/s = 7TeV, 
the energy of the current phase of data taking, and ^/s = 14 TeV, the design energy that will 
be reached in the second phase of operation of the LHC Our selection strategy and the set 
of kinematic cuts that we will design are largely independent, however, of the value of the 
collider energy. This is because they will be optimized to exploit the peculiar kinematics of 
the signal, and a change in the collider energy mainly implies a rescaling of the production 
cross sections of signal and background via the parton luminosities, without affecting the 
kinematic distributions. 

There are two signal topologies contributing to the Wtb channel (see Fig. [s]): either the 
G* decays to one SM and one heavy top (T), and the latter then decays to Wb; or the G* 
decays to one SM and one heavy bottom {B) and the latter decays to Wt. Final states with 
one lepton thus follow from the decay chains: 

qq-^ G* ^ fin + c.c. W^bt + c.c. l^v bb qq 

^ - (16) 
qq-^ G* ^ BbL + c.c. -^W tb + c.c. -^l^vbbqq. 

The physical, observed process is pp — )■ + iets+ ^t, with up to four jets from the hard 
scattering, two of which are 6-jets. The same final state can however follow from another 
signal topology, where G* decays to ti: 

qq^G* ^tt^ W+bW'b ^l^vbbqq. (17) 

This is the decay chain on which previous studies aimed at the discovery of the G* focussed. 

As a reference model for the signal we will consider the two-site model TS5 described in 
the previous section. As explained there, since the decay rates of T — )■ Wbji and B — )■ Wti 



are strongly suppressed in the TS5 (see eq.(14)), the decay chains qq ^ G* ^ Tt+Bb+c.c. — )■ 



Wtb give a negligible contribution and will not be included in the following. We perform 
our analysis adopting the following benchmark values of parameters {nih is the Higgs mass) 

y* = 9*3 = <i , siny^t^ = 0.6 (which implies: tan6'3 = 0.44, sin 2, = 0.55) 

Mg* (18) 
= 1.5 , Mg* = 1.5, 2.0, 3.0, 4.0 TeV , rrih = 120 GeV , 

nif 
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a{pp G*) 
^=7TeV v^=14TeV 


1.5 


1.90pb 


17.4 pb 


2.0 


234 fb 


4.23 pb 


3.0 


4.12fb 


397 fb 


4.0 


0.05 fb 


48.3 fb 



Table 1: Values of the G* production cross section at the LHC for the benchmark choice of 



parameters of eq. ( 18 ) . 



and we make a simple extrapolation to different values at the end. In particular, we fix the 
ratio of the G* mass over the heavy fermion mass to MG*/mf = 1.5, within the interval 
< Mg* < 2mTg/3 where the processes of eq.(16) are kinematically allowed and fast. For 



the reader's convenience, we report in Table [T] the values of the G* production cross section 
at the benchmark points of eq.(18) (see also Fig. [T]). The particles' total decay widths are 



T(G*) T(f) T(B) 

-^-^= 0.052, ^^= 0.035, -^-^ 
Mg* TUf niB 

and the branching ratios of the Wtb and ti channels are |^ 

BR{G* ^ Bb^ Wtb) = 0.17 

BR{G* ^ft^ Wtb) = 0.08 



0.022, 



(19) 



BR{G* ^ti) = 0.06. 



(20) 



Hence, approximately two thirds (one third) of the Wtb events that follow from a heavy-light 
decay of G* is from G* Bb {G* ft). 



The benchmark values of eq.(18) have been chosen so as to have a large rate in the 



heavy-light channels. This particular choice however also implies a large branching ratio to 
a pair of light SM quarks, BR(G* qq) = 0.22 (see Fig. [3]), so that the dijet searches will 
also have a large discovery (or exclusion) power. |^ As previously explained in section 2.2 



this is not generic to the model, and it does not occur for other choices of the parameters, 
like that of the left plot of Fig. |4| which implies large branching ratios for the heavy-light 
channels and a very small one for qq. 



^Here and in the following we use the short notation G* Bb + Tt to denote G* — Bb + ft + c.c 
^In fact, the benchmark point (18 1 might be already excluded, although just marginally, by the bounds 



from the dijet searches reported by the CMS and ATLAS collaborations T2p3 during the completion of this 



T2IT3 



work. This is what suggests a naive rescaling of the bounds on axigluons and colorons shown in Refs 
A more accurate derivation is however required to precisely determine the impact of the dijet searches on 
the parameter space of our model. We leave this to a future study. 
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3.1 Montecarlo simulation of signal and background 



We simulate the signal by using MadGraph v4 |57] , while for the background (see discussion 
below) we make use of both MadGraph and ALPGEN |58] . In our parton- level analysis 
jets are identified with the quarks and gluons from the hard scattering. If two quarks or 
gluons are closer than the separation AR = 0.4, they are merged into a single jet whose 
four-momentum is the vectorial sum of the original momenta. We require that the jets and 
the leptons satisfy the following set of acceptance and isolation cuts: 

PTi > 30 GeV < 5 A% > 0.4 

(21) 

Pti > 20 GeV \r]i\ < 2.5 ARji > 0.4 . 

Here pxj {pti) and rjj (rji) are respectively the jet (lepton) transverse momentum and pseu- 
dorapidity, and ARjj, ARji denote the jet-jet and jet-lepton separations. Jets that fail the 



Pt cut of eq.(21), as well as leptons that fail the isolation and px cut, are discarded and not 
considered in the reconstruction of the event. 6-jets are assumed to be correctly tagged with 
an efficiency ei, = 0.6 if their pseudorapidity satisfies \rih\ < 2.5. We set the corresponding 
rejection rate on light jets to (b = 100. 

Detector effects are roughly accounted for by performing a simple Gaussian smearing on 
the jet energy and momentum absolute value with AE/E = 100%/ E/GeV, and on the jet 
momentum direction using an angle resolution A(f) = 0.05 radians and Arj = 0.04. Moreover, 
the missing energy of each event has been computed by including a Gaussian resolution 
a{^T) = « ■ a/XIj Eij,/GeV , where ^ • Elj. is the scalar sum of the transverse energies of all 
the reconstructed objects (electrons, muons and jets). We choose a = 0.49. 

Table [2] shows the fraction of G* — )■ Bb + Tt events where the jet content is reconstructed 
to be respectively 2j + 26, Ij + 26, Ij + 16 (where j denotes a light jet) for different values 
of the G* mass and ^/s = 14 TeV. Very similar numbers hold for ^/s = 7 TeV. In the third 
column, the numbers in parenthesis show the relative fraction of Ij + 26 events in which 
the light jet originates from the merging of the qq pair from the hadronic W. As naively 
expected due to the increasing boost of all the decay products, the number of events with 
four jets decreases for larger G* masses, while the fraction of events in which the hadronic W 
is reconstructed as a single jet increases. For this reason in our analysis we will select events 
with at least three jets and exactly one lepton and two b-tags: 

pp ^ l^+njets + , n > 3 , 2 6-tags (22) 



where all objects must satisfy the acceptance and isolation cuts of eq.(21). Most likely, the 
use of boosted jet techniques can lead to a better sensitivity on events with three jets, 

-"^"The choice of MadGraph to simulate the signal is dictated by the possibility of easily implementing our 
model. For the background, we have used ALPGEN to simulate processes with a large number of particles 
from the hard scattering. The factorization and renormalization scales have been set to be equal and chosen 
as follows: Q = Mq. for the signal; Q — ^Jm^ + J2j Prj WWbb; Q — yj + p|,^ for Wbbj, Wbbjj, 
WbbSj, WSj and W4j. 

^^The rejection rate corresponds to the inverse of the probability of mistagging a light jet as a 6-jet. 

^^The numerical values of the above parameters, as well as the 6-tagging efficiency and rejection rate, have 
been chosen according to the performance of the ATLAS detector 59 . 

^^See for example the report prepared for the BOOST2010 workshop '60^ for a review and a complete list 
of references. 
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IVIG* 


9t _l_ 97i 


ij + ZO 


ij -|- iO 


1.5 TeV 


0.42 


0.31 (58%) 


0.07 


2.0 TeV 


0.29 


0.42 (79%) 


0.10 


3.0 TeV 


0.13 


0.52 (89%) 


0.17 


4.0 TeV 


0.07 


0.53 (93%) 


0.25 



Table 2: Fraction of G* — t- Bb + Tt events where the jet content is reconstructed to be respectively 
2j + 26, Ij + 26, Ij + 16 (where j denotes a Hght jet) as a function of the G* mass for ^/s = 14 TeV. 
In the third column, the numbers in parenthesis show the relative fraction of Ij + 26 events in 
which the light jet originates from the merging of the qq pair from the hadronic W. 



0.25 




Figure 6: Distribution of the boost factor 7 = E/m of T (left plot) and of the top quark (right 
plot) in signal events G* — )■ Tt with Mq* = 1.5 TeV (continuous orange line) and Mq* = 3.0 TeV 
(dashed blue line) for y/s = 7 TeV. Notice that having set MG*/mj, = 1.5 implies that T is less 
boosted for Mq* = 3.0 TeV than for Mq* = 1.5 TeV. All the curves have been normalized to 
unit area. 



and allow the study of the extreme case, most probable at high G* masses, in which the 
hadronically decayed top is reconstructed as a single fat jet. Notice, on the other hand, 
that while the SM (top or bottom) quark originating from the G* decay is always highly 
boosted, for mG*/mf = 1.5 the heavy quark is not. This is clearly shown by Fig. [6] in the 
case of G* — )■ Tt events. As a consequence, an event selection strategy which relies on a 
large boost of all the decay products, adopted for example by some of the LHC searches for 
heavy resonances decaying to tt, might have a poor efficiency on our signal. 

The largest SM background after the event selection of eq.(22) is the irreducible back- 
ground WWbb, which includes the resonant sub-processes Wtb — )■ WWbb (single top) and 
ti — )■ WWbb. The latter, in particular, gives the largest contribution. We have simulated the 
WWbb events by using MadGraph. Another background which will turn out to be impor- 
tant after imposing our full set of kinematic cuts is Wbb + jets. We have simulated Wbbj, 
Wbbjj and WbbSj (this latter process only for ^/s = 14 TeV) using ALPGEN. Including all 
these samples with increasing multiplicity of light jets in the final state is redundant, and in 
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principle leads to a double counting of kinematic configurations. A correct procedure would 
be resumming soft and collinear emissions by means of a parton shower, and adopting some 
matching technique to avoid double counting. For simplicity, in our analysis we retain all 
the Wbb + njets samples; in this way we expect to obtain a conservative estimate of the 
background. Notice also that some of the cuts we will impose tend to suppress the events 
with larger number of jets and thus to reduce the amount of double counting. Finally, the 
last background that we have considered is + jets, where two light jets are mistagged as 
6-jets. We have generated the processes W3j and WAj with ALPGEN. As for the Wbb + jets 
background, we conservatively include both these sample. In this case, however, the issue 
of double counting can be safely ignored since the W + jets background will turn out to be 
much smaller than the others at the end of our analysis. 

We did not include other reducible backgrounds which are expected to be subdominant, 
in particular: bb + jets where one light jets is misreconstructed as a lepton (it should be 
possible to reduce it down to a negligible level by requiring enough missing energy in the 
event); single-top processes t + jets, tb + jets, Wt + jets (after the request of two 6-tags 
all these are expected to be much smaller than the single-top background Wtb, which is 
included in our analysis). 

While our analysis is carried out at the parton level, we expect it to be robust against 
the inclusion of detector and showering effects. Our simple Gaussian smearing of the jets' 
energy and momentum should correctly reproduce the main impact of detector effects on our 
event selection and reconstruction. The request of two 6-tags and the kinematic cuts that 
we will impose (like for example on the invariant mass of the hadronic W, see next section) 
should reduce the effect of extra radiation in the signal. A detailed study of initial and final 
state radiation, underlying event and multiple parton collisions is beyond the scope of the 
present paper and is left for future analyses. 



3.2 Event selection 



The second column of Tables |3] and |4] reports the value of the cross section for the signal 
and the main SM backgrounds after the selection (22) based on the acceptance and isolation 
cuts of eq.(21), respectively for y/s = 7TeV and 14TeV. One can see that at this stage the 



background dominates by far over the signal. We can however exploit the peculiar kinematics 
of the signal to perform a first set of cuts and reduce the background to a much smaller level 
without touching the signal. 

One basic feature of the signal is that the lepton and the jets in the event tend to be 
very energetic, as they are the final products of the decay chain of a new heavy particle 
(either the G* or the heavy fermion). This is evident for example from the plots of Fig. [7| 
which show the pt of the leading jet (ji) and of the next-to-leading jet (^2) for the signal 
at Mg* = 1.5 TeV and for the total background. We will impose a cut on pT(ji,2) in the 



following, see eqs.(23 ),(24). 



In order to fully exploit the specific topology of the signal, however, it is convenient to 
reconstruct the intermediate Wtb state that follows from the decay of the G* and of the 
heavy fermion, see Fig. [5) To do this, we first reconstruct the longitudinal momentum of 



^See for example Table 32 at page 32 of Ref. 
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LHC Vi 


= 7 TeV 


acceptance 


u + top rec. 


zero-cost 




= 1.5 TeV 


G* ^ft + Bb 
G* ti 


29.8 
5.85 


23.8 
4.72 


22.2 
3.21 


Mr<* 


= 2.0 TeV 


G* ^ft + Bb 
G* tt 


3.29 
0.71 


2.59 
0.56 


2.51 
0.34 




= 3.0 TeV 


G* ^ft + Bb 
G* ti 


0.04 
0.06 


0.03 
0.05 


0.03 
0.02 






WWbb 


4838 


3932 


167 






Wbbj 


210 


156 


7.78 






Wbbjj 


102 


67.5 


2.17 






W3j 


18.8 


14.5 


0.83 






WAj 


8.89 


6.36 


0.31 






Total 












background 


5177 


4177 


179 



Table 3: Cross sections, in fb, at ^/s = 7 TeV for the signal and the main backgrounds after the 



selection (22) based on the acceptance cuts of eq.(21) (second column); after the neutrino and top 



quark reconstruction (third column); after imposing the 'zero-cost' cuts of eq.(23) (fourth column). 
For each channel, the proper branching ratio to a one-lepton final state has been included. 



0.20 
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Figure 7: Differential distribution of the pT of the leading jet (left plot) and next-to-leading jet 
(right plot) after the acceptance cuts of eq.(21) for ^/s = 7 TeV. The continuous line and the dashed 
line respectively show the signal at rn-c* = 1.5 TeV and the total background. All the curves have 
been normalized to unit area. 
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LHC 


= 14 TeV 


acceptance 


u + top rec. 


zero-cost 




= 1.5 TeV 


G* ^ft + Bh 
G — > tt 


262 
42.6 


209 
34.8 


192 
23.7 




= 2.0 TeV 


G* ^ft + Bh 
G* — > tt 


56.9 
7.00 


45.0 
5.76 


42.7 

3.93 


Mg^ 


= 3.0 TeV 


w r J- i ^ uu 

G* tt 


o.ou 
n 46 


38 


9 Q1 
0.22 


Mg* 


= 4.0 TeV 


(_T — f It -t DO 

G* ti 


U.oz 
0.08 


U.Z4 

0.07 


0.03 






WWbb 


27671 


22383 


724 






Wbbj 


794 


573 


25.9 






Wbbjj 


574 


354 


9.30 






WbbSj 


215 


119 


0.63 






W3j 


67.6 


51.3 


2.90 






WAj 


41.2 


28.2 


1.15 






Total 












background 


29363 


23509 


764 



Table 4: Cross sections, in fb, at ^/s = 14 TeV for the signal and the main backgrounds after the 



selection (22) based on the acceptance cuts of eq.(21) (second column); after the neutrino and top 



quark reconstruction (third column); after imposing the 'zero-cost' cuts of eq.(24) (fourth column). 



For each channel, the proper branching ratio to a one-lepton final state has been included. 
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the neutrino, up to a twofold ambiguity, by enforcing the on-shell condition m(/z/) = mw- 
Events where the second-order equation has no (real) solution are removed. We thus have 
two pairs (Zi/), one for each neutrino solution, which represent our two candidates for the 
first W in the signal: PV/i 2 — {^^1,2)- Next, we use the fact that in the signal there is a 
second W that decays hadronically, and label all the jets other than the two 6-jets to form 
our hadronic W candidate, Wh- Notice that in events with exactly three hard jets, Wh will 
consist of a single jet. Starting from the three W candidates and the two 6-jets in the event, 
there are six (Wb) pairs that one can form. We select the pair whose invariant mass is closest 
to the top quark mass, and label it as our top quark candidate: (Wtbt). The remaining b and 
W candidates will be labeled as bf and W^_^ ^ , since in the signal they do not come from the 
decay of a top quark. In the case in which the W selected as belonging to the top quark is 
one of the two leptonic candidates, Wt = Wi., the other leptonic candidate is discarded. In 
other words, the top reconstruction in this case gives us a criterion to select one of the two 
neutrino solutions. If instead the W selected as belonging to the top quark is the hadronic 
one, Wt = Wh, then both the remaining W candidates are kept. As a first loose cut on 
the invariant mass of the reconstructed top we require 80 GeV < m{Wtbt) < 250 GeV. The 
efficiency of this cut on the signal and on the backgrounds with at least one top quark is 
~ 100%. The value of the cross sections after both the neutrino and top reconstruction are 
reported in the third column of Tables |3] and |4j 

The identification of the Wtb intermediate state allows us to fully take advantage of 
the peculiarity of the signal, where each of these three SM particles is extremely energetic. 
This is illustrated by Fig. |8j which shows the distributions of the transverse momenta of the 
reconstructed top, priWibt), and of the W candidate and 6-jet not belonging to it, priW^), 
Pribf), at y/s = 7TeV. Also shown is the distribution of the invariant mass of the hadronic 
W candidate, m{Wh). Similar plots are obtained for y/s = 14TeV. As a first set of cuts, we 
thus use the kinematic observables of Figs. [7| [s] At ^/s = 7TeV we require: 

pAji) > 155 GeV j9t(j2) > 75 GeV m{Wh) < 200 GeV 

(23) 

pAWA) > 105 GeV PriWf) > 90 GeV pt(M > 65 GeV , 

while at ^/s = 14 TeV we impose slightly stronger cuts as follows: 

PtUi) > 175 GeV pr(j2) > 85 GeV m{Wh) < 200 GeV 

PriWA) > 110 GeV Pt{W^) > HO GeV pt(M > ^0 GeV . 



(24) 



In each case, the numerical values of the cuts have been chosen so that each cut individually 
reduces the signal G* Tt + Bb at Mq* = 1.5 TeV by no more than ~ 2%, so that 
they are basically at 'zero cost' for the signal. Notice that when two W^ candidates exist, 
by construction they have the same transverse momentum, so that the cut on priW^) of 



eqs.(23),(24) applies to either of them. As clearly illustrated by Fig. [8| the cut on m{Wh) 
is useful to reduce the backgrounds that do not have a second W , in particular Wbb + jets 
and W + jets. 



-"^^Although this algorithm is very rough and does not account for the amount of off-shellness of the W 
boson, it has a sufficiently large efficiency on our signal: ~ 0.8. We do not attempt here to adopt more 
refined strategies, which are however well known and can be easily implemented in a full analysis. 
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Figure 8: Differential distributions after tlie neutrino and top reconstruction for y/s = 7TeV. 

Upper left plot: pT of the top candidate, priWibt); Upper right plot: px of the W candidate 
labeled as not belonging to the top quark, pTiWf) (if two such candidates exist they have the same 
transverse momentum by construction, see text); Lower left plot: pT of the 6-jet labeled as not 
belonging to the top quark, pT{bf); Lower right plot: invariant mass of the hadronic W candidate, 
m{Wh)- In this latter plot, the first bin is populated by events with exactly three jets, for which 
m{Wh) = in our partonic analysis. The continuous line shows the signal at niG* = 1.5 TeV; in the 
upper plots and the lower left plot, the dashed line shows the total background; in the lower right 
plot, the dashed line shows the (sum of the) backgrounds without a second W, namely Wbb + jets 
and W + jets. All the curves have been normalized to unit area. 
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Figure 9 

mtot = m{WtbtWfbf) , after the cuts at 'zero cost' of eq.(23) at ^/s 



Differential cross section as a function of the total invariant mass of the Wtb system, 

7TeV. The dashed curve 



denotes the total background; the orange, purple, and blue continuous curves denote the signal 
respectively at Mq* = 1.5, 2.0, 3.0 TeV. The long tail of the signal curves at low mtot is due to 
events with an off-shell G* . When two W/ candidates exist in an event, the corresponding two 
solutions for mtot have been both included, each with weight 1/2. 



The cross sections of the signal and the backgrounds after the cuts at 'zero cost' of 
eq.(23) and (24) are reported in the fourth column of respectively Table |3] and |4| At this 
stage, although the background has been strongly reduced, it is still dominant over the signal. 
In particular, the largest background comes from the resonant contribution tt — )■ WWbb. An 
obvious strategy to suppress it, adopted in previous analyses aimed at uncovering the signal 
G* — tt, is that of cutting on the total invariant mass of the event. This is illustrated by 
Fig. [9} which shows the differential cross sections of the signal and of the total background 
for ^/s = 7 TeV as functions of the invariant mass of the Wtb system after the cuts of eq.(23 ). 



Obtaining a better significance of the G* — )■ ti signal over the background, especially at large 
G* masses and widths, requires additional and more sophisticated tools, like for the example 
the use of a left-right polarization asymmetry [5||6]. The peculiar topology of the process 
G* Tt + Bb, on the other hand, suggests a further simple strategy, namely requiring that 
the invariant mass of the (W^bf) system be much bigger than the top mass m^. In the case 
of G* — )■ Tt, indeed, m(Wfy) peaks at mf, and even in G* — )■ Bb events it tends to be much 



larger than m^. This is shown by the contour plot of Fig. 10, which reports the isocurves 
of the (doubly differential) cross section in the plane {mtot,^wb), where mwb = m{Wfy) 
and nitot = miWfyWtbt) . By simple inspection, one can see that the SM ti background 
can be strongly reduced by cutting at the same time on mtot and on mwf,. We find that an 
additional cut on mwt = rniW^Wtbt) is useful at large Mq* to further reduce the background 



and increase the signal significance, see Fig. 11, Although cutting on niwb removes almost 



entirely the tt component of the signal, this can still be detected by adopting the strategies 
proposed in previous analyses. In the following, instead, we will focus on the G* — Tt + Bb 
signal. 
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Figure 10: Contour plot of the signal and total background cross sections after the cuts at 'zero 



cost' of eq.(23) in the plane {mtot,rawb)i where mwb = m{W^h^) and rritot = rn{W^b^Wtbt) ■ The 
center-of-mass energy is y/s = 7 TeV. Different colors correspond to areas of different values of the 
doubly differential cross section (fi a / dmy/bdrntot: as reported in the vertical color key on the right 
in units ab/GeV^. The signals at Mq* = 2.0 TeV and 3.0 TeV have been rescaled respectively by 
a factor 10 and 100. 



For each value of Mq* we find a set of optimized cuts that minimizes the integrated 
luminosity needed for a 5cr discovery. These are: 

Mg* = 1.5 TeV : rutot > 1300 GeV mwb > 400 GeV 
Mg. = 2.0 TeV : rritot > 1700 GeV mwb > 600 GeV 

(25) 

Mg* = 3.0 TeV : mtot > 2500 GeV mwb > 600 GeV mwt > 700 GeV 



Mg* = 4.0 TeV : mtot > 3200 GeV mwb > 700 GeV mwt > 900 GeV . 

In each case, the efficiency of the cut on m-wi, after imposing that on mtot is of the order of 
a few percent for the WWbb background, and ~ 90% for the G* ^ Tt + Bh signal. The 
values of the cross sections after these optimized cuts are reported in Tables |5] and [6] (in 
the columns labeled as 'opt') respectively for ^/s = 7 TeV and 14 TeV. The values of the 
corresponding discovery luminosity are reported in Table [7j 

We define the discovery luminosity to be the integrated luminosity for which a goodness-of-fit test of 
the SM-only hypothesis with Poisson distribution gives a p-value = 2.85 x 10~^, which corresponds to a 5(j 
significance in the limit of a gaussian distribution. If however this value is less than the luminosity for which 
the total (signal plus background) number of expected events is equal to 10, this latter value is conservatively 
defined as the discovery luminosity. 
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LHC 7TeV 


Mg* = 1.5 TeV 

OPT OPT(ii) 


Mg* = 2.0 TeV 

OPT OPT(ii) 


Mg* = 3.0 TeV 

OPT 


G* ^ft + Bb 

WWbb 

Wbbj 

Wbbjj 

WZj 

WAj 

Total 

background 


20.0(1) 14.9(1) 

0.17(6) 0.06(4) 
0.32(2) 0.13(1) 
0.11(2) 0.06(1) 
0.082(3) 0.036(2) 
0.039(3) 0.015(2) 

0.72(7) 0.29(4) 


2.28(1) 1.85(1) 

0.06(4) 0.02(3) 
0.044(7) 0.022(5) 
0.023(8) 0.008(5) 
0.018(2) 0.008(1) 
0.011(1) 0.005(1) 

0.15(4) 0.06(3) 


0.0301(2) 

< 0.03 

0.003(2) 

0.003(4) 

0.0009(3) 

0.0004(4) 

007+"-*^^ 



Table 5: Cross sections, in fb, at ^/s = 7 TeV for the signal and the main backgrounds after 
imposing the cuts of eqs.(21 ),(23), and the optimized cuts of eq.(25) (columns labeled as 'opt'). In 
the case Mg* = 1.5 TeV and 2 TeV the columns labeled as 'opt(ii)' report the value of the cross 



section after the alternative set of optimized cuts of eq.(26), in addition to those of eqs.(21),(23). 
For each channel, the proper branching fraction to a one-lepton final state has been included. 
Details on how the statistical errors and upper limits on the cross sections have been computed 
and combined are given in Appendix IB] 



LHC 14 TeV 


Mg* = 1.5 TeV 


Mg' = 2.0 TeV 


Mg* = 3.0 TeV 


Mg* = 4.0 TeV 


opt 


opt(ii) 


OPT 


opt(ii) 


OPT 


OPT 


G* ^ft + Bb 


175.7(8) 


133.0(7) 


39.5(2) 


32.3(2) 


2.76(1) 


0.231(1) 


WWbb 


3.9(3) 


0.9(2) 


1.0(2) 


0.13(7) 


0.02(5) 


< 0.04 


Wbbj 


2.8(1) 


1.34(8) 


0.76(6) 


0.37(4) 


0.06(2) 


0.005(7) 


Wbbjj 


1.20(4) 


0.59(3) 


0.32(2) 


0.16(1) 


0.028(6) 


0.002(2) 


Wbbjjj 


0.10(1) 


0.037(8) 


0.022(6) 


0.012(5) 


< 0.002 


< 0.002 


W3j 


0.60(2) 


0.26(1) 


0.21(1) 


0.095(8) 


0.043(5) 


0.009(2) 


WAj 


0.28(2) 


0.16(1) 


0.10(1) 


0.055(7) 


0.025(5) 


0.003(1) 


Total 














background 


8.9(3) 


3.3(2) 


2.4(2) 


0.82(8) 


0.18(5) 


(0.0191°:°^,) 



Table 6: Cross sections, in fb, at y/s = 14 TeV for the signal and the main backgrounds after 
imposing the cuts of eqs.(21 ),(24), and the optimized cuts of eq.(25) (columns labeled as 'opt'). In 
the case Mg* = 1.5 TeV and 2 TeV the columns labeled as 'opt(ii)' report the value of the cross 



section after the alternative set of optimized cuts of eq.(26), in addition to those of eqs.(21 ),(24). 
For each channel, the proper branching fraction to a one-lepton final state has been included. 
Details on how the statistical errors and upper limits on the cross sections have been computed 
and combined are given in Appendix [B| 
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Figure 11: Contour plot of the signal and total background cross sections after the cuts at 'zero 
cost' of eq.(23) in the plane {mwt,mwb)i where mwb = m{Wfbf) and mwt = rn{WfWtbt)- The 
center-of-mass energy is y/s = 7 TeV. Different colors correspond to areas of different values of the 
doubly differential cross section cPa /dniwbdmtot, as reported in the vertical color key on the right 
in units ab/GeV^. The signals at Mq* = 2.0 TeV and 3.0 TeV have been rescaled respectively by 
a factor 10 and 100. 



In the case of Mq* = 1.5, 2.0 TeV, we found that imposing cuts shghtly stronger than 



those of eq.(25) can lead to a higher signal over background ratio, S/B, (although the 
corresponding discovery luminosity is also higher). Specifically, we find that the set of cuts 
that maximizes S/B is: 

Mg, = 1.5 TeV : rritot > 1300 GeV mwb > 500 GeV mwt > 900 GeV 

(26) 

Mg* = 2.0 TeV : m^t > 1700 GeV mwb > 600 GeV mwt > HOO GeV . 

The values of the cross sections after imposing this alternative set of optimized cuts are 
reported in the columns of Tables [s] and |6] labeled as 'opt(ii)', respectively for ^/s = 7 TeV 
and 14 TeV. 

3.3 Discovery reach on the parameter space 

All the numbers shown in Tables |5} [6] and [7] hold for the TS5 model at the benchmark point 
of eq.([l8|, where BR{G* ft + Bb ^ Wth) = 0.25 and the G* production cross section 
is that of Fig. [Tj It is however legitimate to ask how these results change when varying the 
model's parameters. The production cross section scales with (tan6'3)^, while the branching 
ratio BR{G* Tt + Bb — )■ Wth) is controlled by tan 6*3, simptji, and the ratio of heavy 
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LHC Vs = 7TeV 


1.5 


Mg' [TeVJ 
2.0 3.0 




Ldisc [fb 


0.48 


4.1 1.3 -lO^ 




LHC Vs= 14TeV 




MG.[TeV] 






1.5 


2.0 3.0 


4.0 


Ldisc [fb 


0.054 


0.24 3.4 


57 



Table 7: Value of the integrated luminosity required for a 5o" discovery after the optimized cuts of 
eq.(25), for = 7TeV (upper panel) and ^/s = 14TeV (lower panel). The 5cr discovery luminosity 



has been computed as explained in footnote 16 



masses Mq* /mf. To simplify the picture, we can fix the latter two parameters to the values 
adopted in the analysis, Mq* /mf = 1.5 and = 3, and study the dependence on tan ^3 and 
siniptR. Figure 12 shows how the branching ratio varies in the plane (sin y^^K, tan ^3). It is 
possible to estimate how the LHC discovery reach varies with tan6'3 and simptR by simply 
rescaling the numbers in Tables [5} |6] to take into account the change in the production cross 
section and in the branching ratio to the final state Wtb. The result is reported in Fig. [T3j 
The two plots show the region in the plane (MG*,tan6'3) where a 5a discovery is possible 
for the LHC at a/s = 7TeV with L = 10fb~^ (upper plot), and at \/s = 14TeV with 
L = 100 fb~^ (lower plot). Three different degrees of compositeness of the right-handed top 
quark are considered: simptR = 0.6,0.8, 1. We expect that this simple rescaling reproduces 
reasonably well the actual reach one would obtain by optimizing the cuts at each different 
point of the parameter space. A possible bias can arise in the limit in which the decay width 
of the G* or those of the heavy fermions become large, and tight cuts on the invariant masses 
have been applied. In this case, the efficiency of the kinematic cuts will in general depend 
non-trivially on the parameters via the particles' decay widths, and will not be reproduced 
by a simple rescaling. In our analysis the cuts on mwb and mwt are always much below the 
heavy fermions' masses, and also the cut on the total invariant mass should be sufficiently 
below the G* mass to neglect, in first approximation, any deviation from the case of a simple 
rescaling. 

The final results of our analysis, summarized by Tables [5| [6] and [7| and by the plots of 
Fig. 13, are very encouraging. For example, if simftR = 0.6 (a value which almost maximizes 
the branching ratio to Wtb in our model, see Fig. 12), the LHC at ^/s = 7TeV and with 
10fb~^ of integrated luminosity should be able to discover a G* with mass in the range 
Mg* = (1.8 — 2.2) TeV for tan ^3 = 0.2 — 0.5. On the other hand, by running at the 
design c.o.m. energy ^/s = 14 TeV, the LHC discovery reach extends to the mass range 
Mg* = (3.3 - 4.4) TeV for tan 613 = 0.2 - 0.5 with an integrated luminosity L = 100 fb^^ 
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0.4 0.5 0.6 0.7 0.8 0.9 1.0 
Sin((^,R) 



Figure 12: Isocurves of constant branching ratio BR{G* Tt + Bb ^ Wtb) in the plane 
{siinptfj,tan6^) with Mc/rrij, = 1.5 and = 2>. 



4 Discussion 

In order to better quantify the potentiahty of our analysis, it is useful to compare our results 
with those obtained in previous studies. 

Ref. [5j searches for a KK-gluon in the channel G* — > ti and makes use of a LR polariza- 
tion asymmetry to enhance the signal significance after cutting on the total invariant mass. 
The authors consider the case in which the right-handed top quark is fully composite, and 
adopt the following benchmark values of parameters: tan 6*3 = 0.2, sunptR = 1, siinpL = 0.33, 

= 3. The corresponding value of the couplings of G* to the light quarks, top quark and 
bottom quarks is: gc^qq = gc^hnhn = -0.25(3, gctLtL = 9G*bLhL - 93, 9Gnutu = 5 g^. The 
authors then assume that only decays of G* to pairs of SM quarks are kinematically al- 
lowed, and obtain in this way a branching ratio BR{G* — )■ tf) ~ 0.95. At the end of their 
analysis, they find S/B = 2 (1.6) and S/VB = 11 (4.2) for Mq. = 3 (4) TeV, where S (B) 
is the number of signal (background) events expected at ^/s = 14 TeV with an integrated 
luminosity L = 100 fb~^. These numbers are obtained without making use of the LR polar- 
ization asymmetry. This observable was proposed in Ref. (see also Ref. ^) as a way to 
better distinguish the signal from the background, but its quantitative impact on the signal 
significance was not derived by the authors. One way to compare with the above results 
is to consider the signal yield that we predict for tan ^3 = 0.2 and rescale it by a factor 
0.95/ BR{G* Tt + Bb — )■ Wtb), so that both results are normalized in the same way 
{i.e. to the same production cross section and decay branching ratio). In this way we obtain 
S/B = 9.3 (3.0), S/VB = 44 (7.3) for Mq. = 3 (4) TeV with L = 100 fb"^ at = 14 TeV. 
This shows how powerful a search for G* — ?■ Tt + Bb is, as a result of the strong reduction of 

^^In order to be conservative in our estimate, we have included a la upward fluctuation of the background 
from its central value quoted in Table [61 
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V7 = 7 TeV L = 10 fb" 




1.0 1.5 2.0 

Mo- (TeV) 



V7 = 14 TeV L = 100 fb" 




Figure 13: LHC discovery reach in the plane {Mq*, tan 9^). The blue area shows the region 
where a discovery of the signal pp ^ G* ^ Tt + Bh — t- Wth is possible at 5cj with svcupt^ = 0.6, 
Mq* jmf = 1.5 and Y^, = 3. The reach at sin ipt^ = 0.8 and sin (ptj^ = 1 is shown respectively by the 
dashed red curve and the dotted black curve. Upper plot: LHC at -^/s = 7 TeV with an integrated 
luminosity L = 10 fb"^ Lower plot: LHC at ^ = 14 TeV with L = 100 fb~^ 
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the SM background obtained by cutting on the invariant mass of the heavy fermion, compared 
to the 'standard' ti channel, although it must be noticed that the sensitivity of the latter 
is expected to increase when the information on the polarization of the tt pair is included. 
Perhaps, a more sensible way to compare our results with those of Ref. [s] is by rescaling the 
ti signal yield obtained in the latter to the branching ratio BR{G* — )■ ti) predicted in our 
model at Mc'/mf = 1.5. For example we consider the benchmark values of the right plot of 
Fig. |4] (tan 6*3 = 0.2, sinLftR = 0.8, sim^i = 0.41, F^, = 3), for which the branching ratios to 
Wtb and ti are comparable: BR{G* ^ft + Bb^ Wth) = 0.21, BR{G* ^ tt) ~ 0.41. After 
the appropriate rescaling, we find that with L = 100 fb~^ at ^/s = 14 TeV the analysis of 
Ref. 1) in this case predicts S/B = 0.86 (0.69), S/VB = 4.7 (1.8) for Mg* = 3 (4) TeV, while 
our result reads S/B = 2.0(0.66), S/y/B = 9.7(1.6). On the other hand, in the extreme 
benchmark point of Ref. [s], where tn is fully composite, we have checked that the ti channel 
leads to a much better discovery reach than the one predicted by our analysis based on the 



Wtb channel. An opposite result is instead obtained for the benchmark values of eq.(18), 
for which the branching ratio to ti is very small. 

To summarize, we find that as long as the right-handed top quark is not exactly fully 
composite and the heavy gluon can decay to one or two heavy fermions, the branching ratio 
of the ti channel decreases and the analyses based on it become less powerful. If decays 
to two heavy fermions are kinematically forbidden, on the other hand, a search based on 
the heavy-light decays to Wtb is very much competitive for the discovery of G* and possibly 
stronger than the 'standard' ti searches. In addition, it gives the opportunity to discover both 
the heavy gluon and the new fermions at the same time. In fact, the process G* ^ Tt + Bb 
turns out to be competitive for the discovery of the T and B as well. 

It is useful, for example, to compare with the results obtained in Ref. [MJ, where the 
QCD pair production and the EW single production of the heavy bottom B were studied. 
At the end of their analysis, the authors find S/B = 1A5 and a discovery luminosity Ldisc = 
4.3 fb~^ for Mb = 1 TeV at ^/s = 14 TeV. For the same value of the heavy bottom mass, 
Mb = ITeV, assuming MG*/mB = 1.5, our analysis leads to a much smaller discovery 
luminosity, thanks to the larger rate of production of the signal. Taking into account only 
the signal G* — )■ Bb — )■ Wtb (i.e. rescaling the signal yields of Table [6] by a factor ~ 2/3 



to subtract the contribution from T, see eq.([20j)), we find S /B = 5.1 — 30 and a discovery 
luminosity L^isc = (0.63 - 0.068) fb~^ if tan ^3 varies in the range 0.2 — 0.5. This shows 
that the heavy-light decays of G* can be an important mechanism for the production of the 
top and bottom KKs which should be taken into account, together with pair production and 
EW single production, to optimize the experimental searches. 

Besides representing an additional and powerful way to discover the heavy fermions, 
the analysis of the heavy-light decays of G* also gives the opportunity to extract some 
important properties of the new fermionic sector. A thorough analysis of the various heavy- 



light channels such as the one sketched in section 2^ can indeed shed light on the quantum 
numbers of the heavy fermions and allow a determination of their couplings through the 
measurement of the branching ratios of G*. Finding a large coupling of the top partners to 
G*, for example, may give a first indication of their composite nature. 



Noti ce that our estimate of the discovery luminosity is more conservative than the one performed in 
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Ref. 

Table 5 of Ref. 31 for Mb = ITeV, we estimate Ldisc — 6.7 fb~^, instead of 4.3 fb~ 
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, which is based on S/vB. For example, using the signal and background cross sections reported in 
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So far, the experimental searches for new heavy fermions have mostly concentrated on 
their QCD pair production. This is expected to be the main production mechanism in 
the case of a fourth generation of chiral fermions, but it does not necessarily represent 
the best discovery process in the case of new vector-like fermions. At the same time, the 
experimental searches for a new heavy gluon have been so far performed only in the dijet and 
tt channels. Our study shows that the presence of new heavy fermions has a strong impact 
on the phenomenology of the heavy gluon, and opens up new possibilities of discovery via 
its decays to one SM and one heavy fermion. The tt experimental searches that are based 
on a full reconstruction of a pair of top quarks are in fact blind to these new processes. On 
the other hand, searches where the top quarks are not reconstructed, like for example those 
of Refs. 



15 



17 



are sensitive to our signal G* Tt + Bb — )■ Wtb, but they are clearly less 
powerful than an optimized analysis, like the one proposed in this work, which exploits the 
peculiar kinematics of the signal. Optimizing the future experimental analyses to take into 
account the heavy-light decay topologies studied in this paper will give the opportunity to 
better explore the parameter space of theories beyond the Standard Model and possibly lead 
to the discovery of the sector responsible for the breaking of the electroweak symmetry. 



Note Added 



While completing our work we became aware of the analysis of Ref. 62 , where the heavy- 



light decays of a KK gluon are also studied. The importance of these latter to explain the 



forward-backward top quark asymmetry at the Tevatron was originally discussed in Ref. 63 
as it was pointed out to us by J. Santiago. 
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Appendix 

A Lagrangian in the diagonal basis before EWSB 

We collect here a few results on the TS5 model used in section [2j As explained in the text, 
the TS5 lagrangian ^ can be diagonalized, before EWSB, by performing a field rotation 
from the composite/elementary to the mass-eigenstate basis. More details on such rotation 
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will be given in a forthcoming publication f56\ . The diagonalized lagrangian reads 

^ ^gauge ~^ ^fermion ~l~ ^Higgs i.'^'^^ 

r — — r^^*^ 

'~' gauge fiu^ 

+ I (DfilD^Gl - D,GID,G:) + \mI,G';^ + f G,. [G;., G;] (28) 



Cfermion =qipq + + X i^P " "^x) X 

- gs tan 6*3 G*^qYq + 93 (sin^v^^ cot 6*3 - cos^v^^ tan 63) Gl^tpYi^ 

+ 93 '^^7'"'^'^ G;xi^^ + g, (cos^^, cot ^3 - sin^^, tan ^^3) G*^xi'x 
sm 6^3 cos 03 \ A, ^ / 

+ h.c. + 0(s2) 



(29) 



Cmggs = \D^H\^ - V{H) 

+ cos v^L cos v^t/j QlHTr + cos QulHTr - sin cos (ptR qiHTn 

- sin i^tR QuiHtR - Y^ cos v^l sin v^tj?^ QlHIr + Y^ sin sin vj^r qiHtn 

- (s2 sin v?L + S3 cos v^l) [cos ^ptR Q'lHTr - sin (ftR Q'lHtR^ 
+ n Qij^Ti + Y, QuHTl - QrHTl 

+ K cos V3fei?^ QdhHBR + K cos v^fcR Q'lHBr - K sin yj^/j QdbHhR 

- Y^ sin QlHIr - F*S2 cos qhHBR + F*S2 sin vj^^j qiHhR 

- Y^S3 sin v^feR QiHbR + Y^S3 cos QlHBr 

+ n g^ji/^i + n Qdiii^^i + g^i/Si + h.c. 

(30) 

where g = u,d,c,s, ip = tL,bL,tR,bR, and we have defined Q = (T,B), Q' = {T',B'), 
Q„ = (T5/3,T2/3), Qd = (5_i/3,S_4/3), H = ((/.+ >o), H = la^H* = {(pl-^)- X denotes 



any of the heavy fermions, except in the first term in the third line of eq.(29), where it 



denotes a top or bottom heavy partner, T, B, T, B. Finally, S3 and S4 are defined as follows: 

^3 = X2 ^2 —SmCfL, S4 = —2 (31) 
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The coupling of G* to the fermions are thus the following (we neglect terms of 0{s2))'- 



9G*qq = - 5-3 tan 6*3 , 

gG*i,^ =93 (sinVi/> cot 6^3 - cosV*/- tan 6*3) 

sin (f^ cos (f^ 

9g*x^ =93 . n ^ ' 

sin6'3COs6'3 

9g*xx =93 {cos^^x cot O3 - sin^(^^ tan 63) 



q = u,d,c, s (32) 

tjj = tL,bL,tn,bR (33) 

XiJ = TtL,BbL,ftn,Bbn (34) 

X = any of the heavy fermions , (35) 



where sim^t^ = simphL = sin ipi, and we have defined sinyj^^ = sin ipsj^ = simp^, simpf^ = 
simptji, simp^^ = simphji, and simp^ = for any other heavy fermion x- In particular, 
at the benchmark point of eq.(|l8|) we have: gc^gg = -OAAg^, 9G*tRtR = 0-53 g^, gcHLtL = 
9G*bLbL = 9G*bRbii = 0.40 (73. The decay rates of G* to two fermions are: 



r (G* ^ qq) 
V {G* ^ + ^X) 

r (G* ^ XX) 



03 

6 

03 
12 

as 



Mg* tan^03 , 

Mg* {siT?ip^ cot 9^ — cos^<^^tan6'3)^ , 



sinV^ cosV^ 
6 ""^sin^e^cos^e, I M? 



X 



G* 



1 



1 ml 



(36) 
(37) 

(38) 



03 
12^ 



G* 



(cos^V9^ cot ^3 — sin^v9^ tan 6^3) ^ + cot^6'3 



+ 6 (cos^i/?^ cot^6'3 — sin^(y9j^) 



ML 



1-4 



^x 



(39) 



Eq.(36) reports the decay rate to a pair of each species of light fermions, q = u,d, c, s, and 
the contribution of both chiralities has been included. Similarly, eq.(39) reports the decay 
rate to a pair of any of the heavy fermions, x^ and the contribution of both chiralities has 
been included. As above, ip = tL,bL,tji,bfi in eqs.(37),(38). 

The heavy fermions mostly decay to one longitudinally-polarized vector boson or Higgs 
boson plus one top or bottom quark, due to their large Yukawa coupling. The decay rates 
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into the three possible channels are: 



327r 



1 + 



Ml 



1 + 



ml + 2M^ 
Ml 



1 - 2- 



M^^ 



Ml 



(40) 



A 



64n' 



1 + 



Ml 



1 + 



+ 2Mi 



1 



1-2- 



-4 



(41) 



64:71 ^ 



1 + 



Ml 




Ml 



m: 



Ml Ml) M^ 



(42) 



At 0(r), by using the Equivalence Theorem |61j, the couplings Xwx^ -^^x' can be ex- 



tracted from the coefficients of the Yukawa terms in the diagonalized lagrangian ( |27[ ) 
the case of the top and bottom partners, one has: 
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In 



A 



WT 



0. 



AzT = AfeT = cos ipL sin (ptR 



(43) 



AvKB = y* COS siu (^^/^ 



A 



(44) 



Avyf = ^zf = Kf = y* sin V^L cos LptR 



(45) 



-^VKB - KS2COSV2bK, XzB = Kb ~ K siu (y^fefl COS (/JfoR X 0(r) . 



(46) 



Notice, in particular, that the decay of B to Wt is suppressed by a factor S2, since it 
must proceed via the mixing of the elementary bi to the composite B'. On the other 
hand, B can decay to Zbji, hhn at next-to-leading order in r. Since under our assumptions 
S2/sinv?feij = S2/ sinipL = mp/mt < r < 1, it follows BR{B Zb + hb) > BR{I3 -> Wt). 
We thus neglect the decay B — )• Wt in our analysis. 



B Statistical Errors 

In this appendix we give some details on how the statistical errors that appear in Tables [5] 
and |6] are computed and combined. 

Let 0" = A a be the cross section for a given process, where A is the number of events and 
a is a proportionality factor with no uncertainty. Then, if a Montecarlo simulation of the 
process returns n events, the true value A is estimated following a Bayesian approach with 
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See for example the discussion in footnote 2 of Ref. 



27 



31 



flat prior and posterior Poisson probability p{X\n) = A" exp(— A)/r;,!. The latter has mean 
E[X] = n + 1 and variance ^[A] = n + 1, which implies a standard deviation on the cross 
section 6a = a/\/n + 1. All the statistical errors on individual cross sections appearing in 
Tables [5] and |6] have been computed according to this formula. 

In the case in which a set of kinematic cuts is applied to a sample of no simulated events 
with initial cross section ctq, and no event passes the cuts, an upper bound on the final cross 
section can be derived at a given confidence level a as the value A=k for which 

A. 

dX p{X\0) = a. 

It thus follows A=K = — log(l — a), which implies an upper limit on the final cross section 
a < l.l/(no + 1) (To at 68% CL. All the upper limits that appear in Tables |5] and |6] are 
computed according to this formula and are at 68% CL. 

When summing over the cross section of several processes, the statistical errors are com- 
bined in quadrature: 

5 a tot = y/^(^(Xi)2. 

In the case in which on one of the initial cross sections there is an upper bound, cxj < aj at 
68% CL, an asymmetric error on the total cross section is derived as follows: 

The above definition is straightforwardly extended to the case in which more than one of 
the initial cross sections have an upper bound. 
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